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Abstract 
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A. 

Balas  and  Saltzman  fl]  iderttified  several  classes  of  facet  i  ’ucing  inequalities  for  the  three-index  assign¬ 
ment  polytope,  and  gave  O(n^)  separation  algorithms  for  two  of  them.  We  give  O(n^)  separation  algorithms 
for  these  two  classes  of  facets,  and  also  for  a  third  class.  Since  the  three-index  assignment  problem  has 
variables,  these  algorithms  are  linear-time  and  their  complexity  is  best  possible. 
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1.  Introduction 


Consider  three  disjoint  n-sets  and  the  collection  of  all  weighted  triplets  with  one  element  in  each  n- 
set.  The  three-index  assignment  (or  three-dimensional  matching)  problem  asks  for  a  minimum- weight  set 
of  triplets  that  partitions  the  union  of  the  three  n-sets.  It  can  be  stated  as  the  following  0-1  programming 
problem: 

min^{c,jtariyt  :  i  €  /,  j  €  6  K), 

s.t.  :  j  €  J,  i  € /f }  =  1)  V»  G /, 

=  VjG/,  (11) 

1}, 

where  /,  J  and  K  are  disjoint  sets  with  |7|  =  (7|  =  lA"!  =  n.  Let  A  be  the  coefficient  matrix  of  the  constraint 
set  of  (1.1).  Then  iZ=/UJU/fis  the  row  index  set  of  A.  Let  S  be  the  column  index  set  of  A.  Let  Ga  be 
the  intersection  graph  of  A.  Then  S  is  the  node  set  of  Ga-  Let 

P  =  {i  G  ;  Ax  =  e,  X  >  0}, 

where  e  =  (1, ...,  1)^  G  3^"-  Then 

Pi  =  conv{x  G  {0, 1}”  :  X  G  P} 
is  the  three-index  assignment  polytope  of  order  n. 

Balas  and  Saltzman  [1]  described  the  cliques  of  Ga  as  belonging  to  three  di-^tinct  classes  and  showed 
that  cliques  in  two  of  the  three  classes  induce  inequalities  that  define  facets  of  P/.  Furthermore,  they  gave 
an  0(n*)  procedure  for  finding  a  facet-defining  clique-inequality  violated  by  a  given  noninteger  x  G  P,  or 
showing  that  no  such  inequality  exists.  They  also  described  in  [1]  the  odd  holes  of  Ga  and  identified  two 
classes  of  facets  associated  with  odd  holes  that  are  easy  to  generate.  One  class  has  coefficients  of  0  and  1, 
the  other  class  coefficients  of  0,  1  and  2. 

In  Section  2  of  this  paper,  we  give  two  O(n^)  separation  algorithms  for  the  above  mentioned  two  classes 
of  facet  defining  clique-inequalities.  In  Section  3,  we  give  an  O(n^)  separation  algorithm  for  facet  inducing 
inequalities  attributed  with  lifted  odd  holes  of  length  5,  having  coefficients  equal  to  0  or  1.  Since  the  number 
of  variables  of  the  three-index  assignment  problem  is  n®,  the  complexity  of  these  procedures  is  of  the  lowest 
possible  order.  Our  procedures  can  be  used  to  speed  up  the  performance  of  the  branch  and  bound  algorithm 
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for  the  three-index  assignment  problem  developed  in  [2],  which  uses  facet  inducing  inequalities  as  ingredients 
of  a  Lagrangian  relaxation  to  obtain  strong  lower  bounds.  Our  procedures  can  be  extended  to  more  general 
cases.  In  Section  4,  we  present  an  procedure  for  detecting  a  violated  facet-defining  inequality  of 

Pi,  which  is  associated  with  a.  2h+  1-hole  of  Ga  &nd  belongs  to  a  subclass  of  the  first  odd-hole-associated 
class  of  facets  of  [1]. 

We  use  the  same  notation  as  [1],  Especially,  a’  means  a  column  of  A  associated  with  s  £  S.  We  may 
specify  s  =  {i,j,  jfe}  for  i  E  I,  j  E  J,  k  E  K.  For  a  set  (J  C  fZ,  we  use  Qi,Qj  and  Qk  to  denote  the  parts  of 
Q  in  I,  J  and  K  respectively. 

For  literature  on  the  three-index  assignment  problem,  also  see  [3],  ...,  [8]. 

2.  Facets  Induced  by  Cliques 
For  every  s  £  S,  let 

C(s)  =  {<€5;a*a‘>2}. 

Then  the  node  set  C(s)  induces  a  clique  of  size  3n  —  2  in  Ga  [1,  Proposition  2.2].  Such  a  clique  is  called  a 
clique  of  class  2.  Furthermore,  for  n  >  3,  the  inequality 

x(C(s))  =  t  €  C(s)}  <  1  (2.1) 

defines  a  facet  of  Pi  for  every  s  e  5  [1,  Theorem  3.3].  Given  a  noninteger  x  £  P,  an  O(n^)  procedure  was 
presented  in  [1,  Section  4]  to  detect  whether  any  inequality  induced  by  a  clique  of  class  2  is  violated.  We 
now  present  an  O(n^)  procedure  to  do  this. 

Algorithm  2.1.  Suppose  that  x  is  a  noninteger  point  in  P.  Let  </  be  an  integer  greater  than  or  equal 

to  4. 

Step  1.  Let  d,  =  0  for  all  s  £  S. 

Step  2.  Check  x<  for  all  t  e  5.  If 

.  (2.2) 
vn 

then  set  d,  :=  d,  +x<  for  all  s  £  C(i).  If  d,  >  1,  then  stop;  (2.1)  is  violated  by  x  with  this  s.  Otherwise, 
continue. 

Step  3.  For  s  €  5,  if 

d,  >  - — (2.3) 
1/ 

then  check  whether  the  inequality  (2.1)  associated  with  s  is  violated  and  if  so,  then  stop.  Else  continue.  ■ 
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Theorem  2.2.  Algorithm  2.1  determines  in  O(n^)  steps  whether  a  given  x  6  P  violates  a  facet  defining 
inequality  induced  by  a  clique  of  class  2.  The  value  of  v  which  minimizes  the  complexity  of  the  algorithm  is 
6. 


To  prove  this  theorem  we  need  two  lemmas.  The  first  lemma  is  essentially  Lemma  4.1  of  [1].  We  repeat 
it  here  for  convenience: 

Lemma  2.3.  For  any  z  €  P  and  any  positive  number  X,  the  number  of  components  with  value  >  X  is 

<  n/X.  ■ 

Lemma  4.1  of  [1]  states  this  for  1/A  integer,  but  the  result  holds  for  reals  as  well.  The  proof  in  [1]  is 
somewhat  indirect.  In  fact,  the  sum  of  components  of  z  €  P  is  always  n,  and  all  the  components  of  z  are 
nonnegative.  The  conclusion  thus  follows. 

Lemma  2.4.  For  any  z  6  P  and  any  positive  number  X,  the  number  of  s  S  S  such  \at  z(C'(s))  >  A  is 

<  n(3n-2)/A. 

Proof:  Consider 

J^{z(C(s))  ;  s  6  5}. 

For  each  t  €  5,  the  number  of  distinct  s  such  that  t  e  C'(s)  is  3n  -  2.  Thus, 

^{z(C(s))  :  s  G  5}  =  (3n  -  2)  :  <  €  S}  =  (3n  -  2)n. 

Since  z(C(s))  >  0  for  each  s,  the  conclusion  follows.  ■ 

Proof  of  Theorem  2.2:  We  first  show  that  Algorithm  2.1  works.  By  Proposition  2.2  of  [1],  |C(s)|  =  3n  — 2 
for  every  s  E  S.  Suppose  that  (2.1)  is  violated  for  an  s  6  5.  Then 

d,  :=  ^{z,  :  t  €  C7(s),x,  > 

>  1  -  :  I  G  C{s),  X,  <  — } 

'  i/n 

^  3n-2^i/-3 

—  *  .  — — — 
i/n  1/ 

i.e.,  (2.3)  holds  and  the  violation  is  discovered  in  Step  3.  Therefore,  Algorithm  2.1  determines  whether  a 
given  z  €  P  violates  a  facet  defining  inequality  induced  by  a  clique  of  class  2. 
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We  now  consider  the  complexity  of  the  algorithm.  The  complexity  of  Step  1  is  n^.  According  to  Lemma 
2.3,  there  are  at  most  components  of  x  satisfying  (2.2).  For  each  t,  the  number  of  s  in  C(t)  is  3n  —  2. 
Hence,  the  complexity  of  Step  2  is  i/n^(3n  —  2).  Finally,  if  (2.3)  holds,  then 

z(C(s))  >  d.  =  :  t  €  C(s),  X.  >  ^}  > 

By  Lemma  2.4,  the  number  of  s  such  that  (2.3)  holds  is  not  greater  than  The  number  of  basic 

operations  to  check  (2.1)  is  3n  —  2.  We  also  need  to  use  basic  steps  to  check  (2.3).  Therefore,  the 
complexity  of  Step  3  is  n^+  Hence,  the  overall  complexity  of  Algorithm  2.1  is  O(n^).  The  value  of 

1/  that  minimizes  the  maximum  of  the  two  polynomials  expressing  the  complexity  of  Steps  2  and  3  is  easily 
seen  to  be  6.  ■ 

We  now  consider  facets  induced  by  cliques  of  class  3  [1].  Suppose  that  s,t  g  5  such  that 

a’-a‘  =  0.  (2.4) 


Let  s  =  t  =  ti  =  h  =  (»,,;,, Ir,),  <3  =  and 

C(s,t)=  {S,h,t2,t3}- 

Then  the  node  set  C(s,t)  induces  a  (4-)clique  in  Ga  [1,  Proposition  2.3].  Such  a  clique  is  called  a  clique 
of  class  3.  The  number  of  cliques  of  class  3  is  \n^{n  —  1)®  [1,  Corollary  2.5).  Furthermore,  for  n  >  4,  the 
inequality 

x(C(s,  <))  =  ;  t  e  C{s,  /)}  <  1  (2.5) 

defines  a  facet  of  Pi  for  every  s,(  €  S  satisfying  (2.4)  (1,  Theorem  3.7].  Given  a  noninteger  x  G  P,  an  0(n‘*) 
procedure  was  also  presented  in  [1,  Section  4]  to  detect  whether  any  inequality  induced  by  a  clique  of  class 
3  is  violated.  We  now  present  an  O(n^)  procedure  to  do  this. 


Algorithm  2.5.  Suppose  that  x  is  a  noninteger  point  in  P.  Check  x,  for  s  G  S.  If 


1 

-  <  X,  < 

4 


1. 


then  check  Xp  for  p  G  5  satisfying 

a*  •o'’  =  1. 


If 


Xp  > 


I  -  X, 
3 


(2.6) 


(2.7) 


(2.8) 
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then  for  t  £  5  such  that  a*  •  a*  =  0  and 

a‘  ■  flP  =  2,  (2.9) 

check  whether  (2.5)  is  violated  or  not.  ■ 

Theorem  2.6.  Algorithm  2.5  determines  in  O(n^)  steps  whether  a  given  x  £  P  violates  a  facet  defining 
inequality  induced  by  a  clique  of  class  3. 

Proof:  We  first  show  that  Algorithm  2.5  works.  Let  C{s,t)  be  the  node  set  of  a  clique  of  class  3.  Since 
|C(s,<)l  =  4,  if  a:  G  P  violates  (2.5),  then  x  has  at  least  one  component  >  By  Proposition  2.4  of  [1],  there 
is  no  loss  of  generality  in  assuming  that  this  happens  for  the  components  indexed  by  s,  i.e.,  that  x,  >  ^. 
Since  x  G  P,  we  have 


:jeJ,k£K}  =  l, 

(2.10) 

:ieI,k£K}  =  1, 

(2.11) 

:  »€/,;■  G  =  1. 

(2.12) 

If  X,  =  1,  then  Xp  =  0  for  any  p  such  that  a*  =  1,  by  (2.10),  (2.11)  and  (2.12).  In  this  case,  for  any  t  such 
that  a‘  •  a*  =  0,  the  left  hand  side  of  (2.5)  is  1.  Thus,  if  x  G  P  violates  (2.5),  x,  <  1.  This  justifies  (2.6). 
Since  |C(s,<)|  =  4,  if  (2.5)  is  violated,  then  (2.8)  must  hold  for  at  least  one  p  G  C(s,t),  p  ^  s.  Therefore, 
this  algorithm  determines  whether  a  given  x  G  P  violates  a  facet  defining  inequality  induced  by  a  clique  of 
class  3. 

We  now  consider  the  complexity  issue.  The  complexity  to  check  all  s  is  0{n^).  According  to  Lemma 
2.3,  there  are  at  most  4n  components  of  x  satisfying  (2.6).  For  each  such  component  x,,  there  are  (n  -  1)^ 
components  Xp  satisfying  (2.7).  So  the  complexity  to  find  all  pairs  s,p  G  S  satisfying  (2.6)  and  (2.7)  is 
O(n^).  By  (2.10),  (2.11)  and  (2.12),  for  given  x,  satisfying  (2.6),  there  are  at  most  6  p  such  that  (2.7)  and 
(2.8)  hold.  When  s  and  p  are  fixed  and  (2.7)  is  satisfied,  there  are  n  -  1  t  satisfying  (2.4)  and  (2.9).  One 
needs  3  additions  and  one  comparison  to  check  (2.5).  Therefore,  the  complexity  of  the  remaining  part  of  the 
algorithm  is  O(n^).  Hence,  the  overall  complexity  of  Algorithm  2.5  is  0{n^).  ■ 

In  the  proof,  we  see  that  the  main  work  is  to  check  all  s  £  S.  This  work  cannot  be  reduced.  So  this 
order  cannot  be  further  improved.  Similarly,  the  order  in  Theorem  2.2  is  also  the  minimum. 

By  [1],  only  these  two  classes  of  cliques  induce  facets  of  P/.  Combining  Theorems  2.2  and  2.6,  we  have 

Theorem  2.7.  One  can  determine  in  O(n^)  steps  whether  a  given  x  £  P  violates  a  facet  defining 
inequality  induced  by  a  clique  of  G a-  ■ 
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3.  Facets  Associated  with  Five-Holes 


Assume  that  n  >  3.  Let  Q  C  72,  |Q|  =  2A  +  1  for  some  positive  integer  h  with  /i  <  n  —  1, 1  <  \Ql\  < 


h,  L  =:  I,J,K,  and  let 


Then  the  inequality 


5(g)  1=  {se5:5^{aj  :,€«}>  2}. 


;  s  e  5(Q)}  <  h 


defines  a  facet  of  Pi  [1,  Theorem  6.1].  These  facets  can  be  regarded  as  lifted  from  odd  hole  inequalities.  The 
number  of  distinct  inequalities  (3.2)  is  0(2^")  [1,  Proposition  6.5).  In  particular,  facet  defining  inequalities 

(2.1)  are  special  cases  of  (3.2)  with  h  =  1. 

Denote  the  set  of  Q  satisfying  the  above  requirements  with  a  fixed  h  by  5^ .  Let  =  {5(g)  :  g  6  5*  }, 
where  S{Q)  is  defined  by  (3.1).  Thus,  Algorithm  2.1  is  an  O(n^)  algorithm  to  detect  violated  facets  induced 
by  sets  in  <ti.  In  this  section,  we  present  an  O(n^)  procedure  to  detect  whether  there  exists  a  facet  induced 
by  a  set  in  $2.  which  is  violated  by  a  given  noninteger  x  E  P-  Notice  that  the  number  of  facets  defined  by 

(3.2)  with  /i  =  2  is  0(n®).  Without  loss  of  generality,  we  may  assume  that  such  an  x  does  not  violate  any 
inequality  (2.1),  i.e.,  assume  that  Algorithm  2.1  was  applied  to  x  first. 

Suppose  that  x  E  P.  For  iEl,j€J,kEK,  we  will  denote 

x(i,  j,  K)  :=  ^{Xijk'  :  k'  E  K},  (3.3) 

x{i,J,k)  -feJ},  (3.4) 

x{IJ,k):='^{xi.jk  i'  El]  (3.5) 

In  the  following,  we  establish  some  formulas  to  calculate  x(S(Q))  for  Q  E  S^-  We  consider  two  cases. 


First,  suppose  that  s,t  E  S  satisfy 


a*  o'  =  1. 


g  =  sUt  e  5^. 

Assume  that  s  =  {i,,j$,k,),t  =  (it,jt,kt)-  It  is  not  difficult  to  see  that  if  it,  =  ik,,  then 

x(5(g))  =  x(C(s))  +  x(C(<)) +  *(».,;«.  A)  +  *(»■»>;.>  A)  -  Xi.j.t.  -  x<,;,t,; 

if  =  jt,  then 


x(5(g))  =  x{C(8))  +  x(C{t))  +  x(.i,,J,k,)  +  x(it,J,k,)  -  x,.j,fc,  -  Xi,,.*.; 
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and  if  i,  =  I't ,  then 


x(SiQ))  =  x(C{s))  +  x(C{t))  +  z{l,j,,k,)  +  x{I,j,,k,)  -  Xi,j,k,  -  xi,j,t,. 


(3.8) 


Next,  suppose  that  s,t  £  S  satisfy 

a* -a*  =  2.  (3.9) 

Again  assume  that  s  =  {i,,j,,k,),t  =  (it,jt,kt).  Then  define  for  the  ordered  pair  (s,t), 

x(s,t)  :=  x(i,,J,k,)  +  x{it.jt,  K)  -  (3.10) 

if  »i  7^  *«; 

i(s,t)  :=  x(it,jt,I<)-k  x(ljt,k,)  -  2x„,j,,t,,  (3.11) 

if  js  jt ; 

x(s,t)  ;=  x{I,jt,k,)  +  x(it,J,  k,)  -  2x„j,,i,,  (3.12) 

if  k,  ^  kt-  Notice  that  in  (3.10)  jt  =  j,  and  kt  =  k,,  in  (3.11)  it  =  i,  and  kf  =  k,,  and  in  (3.12)  i|  =  i, 
and  jt  =  jf  Thus,  (3. 10)  can  be  written  with  jt  and  k,  substituted  for  jt  and  kt,  and  there  are  similar 
expressions  for  (3.11)  and  (3.12). 

In  the  case  that  (3.9)  holds,  we  may  expand  s  U  t  to  a  set  Q  €  5*  by  adding  an  adequate  element.  For 
example,  if  i,  ^  it ,  we  may  add  j  ^  j,  and  let 

Qi  =  Qj  =  {j.J],  Qk  =  {k.},  (3.13) 

or  we  may  add  k  ^  k,  and  let 

Qi  =  {«».it},  Qj  =  Qk  =  {k,,k}.  (3.14) 

Notice  that  there  are  2(n  -  1)  ways  to  expand  a  given  pair  of  s  and  t  satisfying  (3.9)  to  a  set  Q  e  S^.  We 
may  also  calculate  x{S{Q))  by  expressions  similar  to  (3.6)-(3.8).  For  example,  in  case  (3.13),  one  can  check 
that 


x{S{Q))  =  x(C(s))  +  x(s,t)  +  x(i,,j, /f)  +  i(it, ;,/(■)  + x(/,j,fc,)  -  2xi.jt,  -  2x„ji.,  (3.15) 

where  x(s,t)  is  given  by  (3.10).  Note  that  if  i,  ^  i|,  there  are  two  possibilities,  namely,  (3.13)  and  (3.14). 
Therefore,  there  are  six  cases  totally,  i.e.,  two  cases  for  i,  ^  it,  two  cases  for  j,  ^  jt,  and  two  cases  for 
k,  ^  kt-  In  each  of  these  six  cases,  there  is  an  expression  of  type  (3.15)  to  calculate  x(S(Q)).  Later,  we  will 
refer  to  these  as  “expressions  of  type  (3.15)”. 
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We  are  now  ready  to  describe  our  algorithm; 


Algorithm  3.1.  Suppose  that  x  is  a  noninteger  point  in  P.  Suppose  (2.1)  holds  for  all  s  £  S. 

Step  1.  Let  d,  =  0  for  all  s  S  S. 

Step  2.  Check  X(  for  all  t  E  S.  If 

>  -i,  (3.16) 

then  let  d,  :=  d,  +  X(  for  all  s  €  5  satisfying 

a,  ■  a«  >  2. 


Step  3.  For  all  s  €  5,  if 

(3.17) 

then  set  d,  :=  x(C(s)). 

Step  4.  For  all  i  E  I,j  €  d,  it  €  K,  calculate  x{i,  j,  K),  x{i,  J,  k)  and  x{l,j,k)  by  (3  3),  (3.4)  and  (3.5). 
Step  5.  For  all  i  E  I,j  €  J,k  E  A',  form  the  sets 

L(i,J.  K)  :=  {s  E  S  :  d,  >  ^,i,  =  »}, 

L{I,j,K):=  {sES:d.>  5.7.=;} 

and 

L(I,J.k):=  [sES  :d.>  ^it,  =  F} 

Store  them. 

Step  6.  For  each  s  =  {i,,j,,k,)  E  S  such  that  d,  >  i,  check  dt  for  each  t  E  A(s),  where 

L{s)  :=(L{i.,J,l<)UL{lJ.J<)UL(I,J,k,))\{s}.  (3.18) 

If  a‘  •  a’  =  1  and  d*  >  (2  —  d,)/3,  calculate  x(5(Q))  by  the  appropriate  expression  (3.6),  (3.7)  or  (3.8), 
with  dp  substituted  for  x(C(p)),p  =  8,t. 

If  a'  ■  a'  =  2  and  x(s,t)  >  (2  —  d,)/3,  calculate  x(5(Q))  by  the  appropriate  expression  of  type  (3.15), 
with  dp  substituted  for  x(C(p)),p  =  s,<,  for  every  Q  E  such  that  sU<  C  Q. 

If  x(5(Q))  >  2  for  some  Q,  stop:  the  corresponding  inequality  (3.2)  is  violated  by  x.  Otherwise, 
continue.  ■ 


8 


Theorem  3.2.  Algorithm  3.1  determines  in  O(n^)  steps  whether  a  given  x  £  P  violates  a  facet  defining 
inequality  (3.2)  with  h  =  2. 

We  prove  several  auxiliary  results. 

Lemma  3.3.  After  Step  3  of  Algorithm  3.1,  for  all  s  E  S, 

d,<x{C{s)).  (3.19) 

Ifx(C{s))  >  i,  then  rf,  =x(C(s)). 

Proof:  It  is  obvious  that  (3.19)  holds.  Suppose  that  z(C(s))  >  5.  Then 

\  ~  ^  ^ 

1  3n  -  2  1 

-  3  12n  -  12’ 

i.e.,  d,  >  -^  holds  and  in  this  case  d,  is  set  equal  to  x(C(<)  in  Step  3.  ■ 

Lemma  3.4.  At  Step  6, 

|L(s)|  <  27n,  (3.20) 

|{«€  L(s):a^-a‘=2,x(s,t)>  <  18-  (3-21) 

where  L(s)  is  defined  by  (3.18)  and  x(s,t)  is  defined  by  (3.10),  (3.11)  and  (3.12). 

Proof:  Suppose  that  s  =  ii,,j3,h,). 

On  the  one  hand,  by  Lemma  3.3  and  the  definition  of  L{i,,  J.  K)  at  Step  5, 

J2{^iC(t))  .t€L(i„J,K)}  >  ^|L(.-,,y,A')(. 

On  the  other  hand, 

^{x(C(<))  ;  t  E  L(i,,J,  A')} 

ii.ki 

ii.ki  i  j  * 

=  1351  +  13  51  +  13  5Z 
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=  5;i  +  5^1  +  ^l  =  3n. 

«  ii  ti 

Thus,  \L{i,,J,K)\  <  9n.  Similarly,  we  can  show  that  \L(I,j,,K)\  <  9n  and  \L(I,j,,K)\  <  9n.  From  the 
definition  of  L(s),  we  then  have  (3.20). 

It  is  also  seen  that 

^{x(s,t)  :t€S,a‘  a*  =  2} 

—  2  ^  ^i,ik  +  2  ^ij,k  +  2  ^ijk,  —  6. 

},k  i.k  ij 

Hence,  (3.21)  holds.  ■ 

Lemma  3.5.  Suppose  that  Q  £  S^.  Let 

H{Q)  ■.=  {s€S  :  9  e  Q}  =  3}. 

If  the  inequality  (3.2)  with  h  =  2  is  violated  by  x,  then  there  are  s,t  £  H{Q),  s  such  that 

(i) .  x(C(s))  >  i- 

(ii) .  either  a’  •  a*  =  1  and  x{C{t))  >  (2  -  x(C(s)))/3,  or  o*  ■  a*  =  2  and  x(s,  1)  >  (2  —  x(C(s)))/3. 

Proof:  Without  loss  of  generality,  assume  that  Qj  =  {ji,»2},  Qj  =  {jiji},  Qk  =  {'t}-  Then  H{Q)  = 
{(iijuk),  (iij2,  k),  (h.Ji.k),  (hJjrk)},  and 

2  <  x(5(g))  <  x(C((ii ,  ji ,k)))  +  x(C((ii ,j2,k)))  +  x(C((i2,ji ,k)))  +  x(C((i2 ,j2,k))). 

Then  at  least  one  term  on  the  right  hand  side  of  the  above  inequality  is  greater  than  j.  Without  loss  of 
generality,  assume  it  is  the  first  term.  Let  s  =  (ii,ji,k).  Then  (i)  holds. 

Denote  =  (i2,j2,k),t2  =  {iij2,k),t3  =  (»2.ii.*)-  Then 

2  <  x{S(Q)) 

<  x{C{s)  +  x(C(ti))  +  x(i:,  j2,  K)  +  x{i2,ji,K)  -  xj,;,*  -  Xi,;,*,  (by  (3.6)) 

<  x(C(s))  +  x(C(ti))  +  x(«,t2)  +  *(«.<3)-  (by  (3.11)  and  (3.10)) 

Notice  that  a*  •  a*'  and  a*  ■  a*’  =  a*  ■  a*’  =  2.  Thus, 

2  —  x(C(s))  <  x(C(ti)  +  x(s,  <2)  +  x{s,  <3). 
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At  least  one  term  at  the  right  hand  side  of  above  inequality  is  greater  than  one  third  of  the  left  hand  side  of 
this  inequality.  Hence,  (ii)  holds.  ■ 

Proof  of  Theorem  3.2:  By  Lemmas  3.5  and  3.3,  if  the  inequality  (3.2)  with  A  =  2  is  violated  by  x  for 
some  Q,  there  are  s,t  £  S  such  that  d,  >  and  either  a'  •  a*  =  1  and  dj  >  (2  —  d,)/3,  or  a'  a*  =  2  and 
x{s,t)  >  (2  -  d,)/3. 

Thus  if  t  €  L(s),  i.e.,  if  dt  >  g,  then  Step  6  will  detect  the  violation.  To  see  that  this  is  the  case,  notice 
that  if  •  a‘  =  1,  then 

d«>(2-d.)/3  =  (2-x(C(s)))/3> 

(since  x(C(s))  <  1  for  all  s  €  S),  and  if  a‘  a*  =  2,  then 

dt  =  x(C(t))  >  x(s,t)  (from  the  definitions) 

>(2-d,)/3  =  (2-x(C(s)))/3>  i. 

Thus  in  both  cases  t  €  L(s). 

Consider  now  the  complexity  of  Algorithm  3.1.  Steps  1  and  2  are  obviously  O(n^).  In  Step  3,  the 
number  of  s  €  5  such  that  d,  >  ^  is  <  12n(3n  —  2)  (From  Lemma  2.4  and  the  fact  that  d,  <  i(C(s)); 
and  for  each  number  s,  calculating  x(C(s))  requires  3n  —  3  additions;  thus  the  complexity  of  Step  3  is  also 
O(n^).  Further,  each  of  Steps  4  and  5  again  requires  0(n®)  operations. 

Finally,  consider  the  complexity  of  Step  6.  By  Lemmas  3.3  and  2.4,  the  number  of  s  such  that  d,  >  L 
is  O(n^).  From  Lemma  3.4,  the  cardinality  of  L{s)  is  0(n).  Forming  L{s)  itself  is  0(1),  since  its  constituent 
sets  have  been  calculated  and  stored  in  Step  5.  Thus,  checking  d(  for  each  t  €  L{s)  for  all  s  such  that  d,  >  | 
is  O(n^).  If  a'  •  a*  =  1  and  d(  >  (2  -  d,)(Z,  calculating  x(5(Q))  is  0(1),  since  the  expressions  x(C(s)), 
x(C(i))  and  x{i,j,K),  x{i,J,k),  x{I,j,k)  have  been  calculated  in  Step  3  and  Step  4  respectively.  On  the 
other  hand,  from  Lemma  3.4  the  number  of  those  t  €  L(s)  (for  given  s)  such  that  a*  •  a'  =  2  and  x(s,  <)  >  g 
is  0(1).  For  each  such  t,  there  are  2(n  —  1)  ways  to  expand  the  pair  {s,t)  to  a  set  Q  G  5^,  and  for  each 
such  set  Q,  calculating  x{S{Q))  is  again  0(1),  since  the  constituent  terms  of  the  expressions  of  type  (3.15) 
needed  for  this  have  already  been  calculated.  Thus,  the  complexity  of  dealing  with  all  pairs  (s,t)  satisfying 
a*  ■  a‘  =  2  is  again  O(n^),  which  leaves  the  total  complexity  of  Step  6  at  O(n^).  ■ 

4.  Facets  Associated  with  Other  Odd  Holes 

It  is  somewhat  more  complicated  to  detect  violated  facets  defined  by  (3.2)  with  h  >  2.  We  first  discuss 
some  subclasses  of  5*  and  Assume  that  n  >  3.  Let  S'*  be  the  subcollection  of  5*  such  that  for  any 


11 


Q  E  s''  there  exists  an  L  =  I,  J,  or  K  such  that  \Qi,\  =  1.  Note  that  in  this  case  |Qi'|  =  h  for  L'  G  {/,  <7,  A'}, 
U  ^  L.  Let  =  {5(Q)  :  Q  G  5*}-  Clearly,  =  5‘,  5^  =  5*,  but  §*'  is  a  proper  subset  of  S'*  for  h  >  2. 
Assume  that  h  >  3.  We  will  show  that 

Theorem  4.1.  For  a  given  x  E  P,  H  is  possible  in  0(n^'*’^)  steps  to  detect  whether  there  exists  a  facet 
of  Pi,  which  is  defined  by  (3.2)  with  Q  E  §''  and  violated  by  x. 

Again,  we  prove  several  auxiliary  results. 

Suppose  that  x  E  P  does  not  violate  any  inequr“* ' 

x(S(Q))  <  r.  (4.1) 

with  Q  E  S'" ,r  =  1,2,  1,  but  violates  (4.1)  with  h  =  r  and  a  set  Q  G  S'*.  Suppose  that  h  is  fixed,  not 

dependent  on  n. 

Lemma  4.2.  Let  H(Q)  be  defined  as  in  Lemma  3.5.  Then  there  exists  an  s  E  H{Q)  such  that 

x(C(s))  >  i.  (4.2) 

Proof:  Clearly, 

h  <  x{sm  <  ^ 

However,  for  Q  G  S'*,  \H{Q)\  <  h^.  The  conclusion  thus  follows.  ■ 

Since  Q  E  S*,  without  loss  of  generality,  assume  that  Qk  —  {ife}.  Then  \Qi{  —  \Qj\  =  h.  Suppose  now 
that  Si  =  (»! ,  ji ,  it)  G  H{Q)  and  (4.2)  holds  for  s  =  si. 

Lemma  4.3.  One  of  the  following  clatms  holds. 

(i).  Q,  =  {ti,«2,  ...,ih},  where  i\,i2,—,ih  are  distinct,  and  Qj  3  {juj2,..jT},  where  r  <  h, 
are  distinct,  such  that 

x(C(s,,))  >  _  j  I  (d-3) 

forv  =  2,...,h,  where  s„  =  (iy,j„,*),  and  E  {ji,  .  ^jr]  for  r  <  j  <  h; 

(»)■  Qi  2  . ‘t},  where  r  <  h,  11,12,. ..,Jt  ore  distinct,  and  Qj  =  {ii,j2,  where 

ji,h<  ■■■ijh  are  distinct,  such  that  (4-3)  holds,  for  1/  =  2, ..., h,  where  s„  =  (iv,ji>tk),  and  i„  E  {ii, ...,  ir)  for 
T  <  i  <  h. 
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Proof:  Suppose  that  (4.3)  holds  for  s„,  i/  =  2, r  such  that  there  are  no  repetitions  in  {r’l , i't}  and 
{ii.  Suppose  that  such  a  collection  is  maximal,  i.e.,  for  any  »  €  Q/\{m,  •••,  *t}  and  j  E.  Qj\{ji,  ...Jt)  , 

*(C((.-.i.ib)))<  2^.  (4.4) 

If  the  conclusion  of  this  lemma  does  not  hold,  then  there  exists  I'o  G  0/  \  {»i ,  »V}  and  jo  €Qj\  { ji ,  jr} 

such  that 

(4.5) 

x(C((i.,jo,fc)))<  2^.  (4.6) 

for  u  —  1, ...,  T.  However, 

x{S{Q))  <  x{S{Q‘))  +  5^{x(C(io,  j,ib) :  j  6  Qj,  j  9^  jo)  +  5]{x(C(i,  jo,  k)  :  j  e  Qj),  (4.7) 

where  Q'  is  defined  by  Q'j  =  Qi\  {io}.  Q'j  =  Qj\  {jo},  Q'k  =  Qk  =  {k}.  By  (4.4),  (4.5)  and  (4.6),  each 
term  in  the  right  hand  side  of  (4.7),  except  the  first  term,  is  not  greater  than  number  of  such 

terms  is  2A  -  1.  Clearly,  Q'  €  5^“'.  Thus,  x(S(Q'))  <  h  —  I  by  our  assumption.  Thus,  the  right  hand  side 
of  (4.7)  is  not  great  than  h.  This  forms  a  contradiction  with  our  assumptions  on  x  and  Q.  Therefore,  this 
lemma  holds.  ■ 

The  following  is  a  brief  proof  of  Theorem  4.1.  Cumbersome  details  are  omitted. 

Proof  of  Theorem  f.l:  By  changing  12  to  4(2/»—  1)  in  (3.16),  and  (3.17),  we  may  calculate  all  x(C(s))  > 
in  O(n^)  steps.  Check  s  €  5  with  x(C(s))  >  j.  The  number  of  such  s  is  O(n^).  Let  si  =  s  =  (»j,  ji,  ii). 
We  first  check  Q  E  such  that  s  CQ  and  Q  has  a  structure  as  described  in  Lemma  4.3  (i). 

Let  k  =  ki.  For  r  =  2,...,h,  find  all  possible  =  (ii,,jv,k)  for  u  =  2,...,r  such  that  (4.3)  holds  and 
there  are  no  repetitions  in  {»j,...,iT}  and  {ji,...,jr}.  (We  may  find  S2,...,St  sequentially.)  Similarly  to 
Lemma  3.4,  one  may  see  that  the  number  of  t  =  k)  such  that 

>  5rrr 

is  0(n).  Thus,  the  number  of  such  (s2,...,Sr)  is  0(n’'“*).  Then  we  find  all  possible  s„  =  {iv,jv,k)  for 
1/  =  r  +  1, ...,  h  such  either  (i)  or  (ii)  of  Lemma  4.3  holds.  With  an  argument  similar  to  Lemma  3.4,  one  may 
see  that  the  number  of  such  s„  is  0(1).  In  order  to  expand  si  U  ...  U  s*  to  a  set  Q  G  5^,  we  need  to  add  it 
with 

{jr+l,...,jA}  C  J\  {jl,...,jr}. 
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Therefore,  the  number  of  expansion  ways  is  For  each  Q  produced  in  this  way,  we  may  use  formulas 

simiiar  to  (3.6)-(3.8),  (3.15)  to  calculate  x(S{Q))  and  check  whether 

x(SiQ))  <  h.  (4.8) 

The  complexity  of  such  formulas  are  0(1).  Summing  up  the  above  discussion,  one  see  that  the  complexity 
of  check  (4.8)  for  all  Q  G  such  that  s  C  Q  and  Q  has  a  structure  as  described  in  Lemma  4.3  (i),  for  a 
fixed  s,  is  0(n'*“^). 

The  case  of  Lemma  4.3  (ii)  and  the  cases,  where  i\  or  ji  in  a  is  fixed,  are  similoj^.  Si^'ce  the  number 
of  such  s  is  O(n^),  the  overall  complexity  of  such  a  procedure  is  0(n*''‘*).  This  completes  the  proof  of  this 
theorem.  ■ 

In  general,  the  order  of  procedures  to  detect  violated  facets  induced  by  sets  in  for  h  >  3  may  be 
higher  than  0(n*‘*'*).  For  example,  let  h  =  A.  We  may  find  0(n®)  triples  of  (s,t,p)  such  that  s,t,p  €  S, 
x(C(s)),x(C(<)),x(C(p))>i,and 

a*  ■  a*  =  0,  a‘  ■  a'’  =:  0,  a‘  ■  a’^  =  0. 

Then  0:=sU<Up€5‘‘.  We  need  to  check  whether 

x(sm  <  4 

or  not.  Thus,  the  complexity  of  such  a  procedure  is  at  least  0(n®).  But  6  >  h  + 1  =  5.  Hence,  the  complexity 
for  detecting  violated  facets  induced  by  sets  in  may  be  higher  than  the  complexity  corresponding  to 
for  A  >  3. 
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Balas  and  Saltzman  (1)  identified  several  classes  of  ^acet  inducing  inequalities 
for  the  three-index  assignment  polytope,  and  gove  O(n^)  separation  algorithms  for 
these  two  classes  of  facets,  and  also  for  a  third  class.  Since  the  three-index 
assignment  problem  has  n^  variables,  these  algorithms  are  linear-time  and  their 
complexity  is  best  possible. 
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